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ABSTRACT 

Aims. The primary objective of this study is to search for and identify wave modes within a sunspot penumbra. 
Methods. Infrared spectropolarimetric time series data are inverted using a model comprising two atmospheric components in each 
spatial pixel. Fourier phase difference analysis is performed on the line-of-sight velocities retrieved from both components to de- 
termine time delays between the velocity signals. In addition, the vertical separation between the signals in the two components is 
calculated from the Stokes velocity response functions. 

Results. The inversion yields two atmospheric components, one permeated by a nearly horizontal magnetic field, the other with a 
less-inclined magnetic field. Time delays between the oscillations in the two components in the frequency range 2.5 - 4.5 mHz are 
combined with speeds of atmospheric wave modes to determine wave travel distances. These are compared to expected path lengths 
obtained from response functions of the observed spectral lines in the different atmospheric components. Fast-mode (i.e., modified 
p-mode) waves exhibit the best agreement with the observations when propagating toward the sunspot at an angle ~50° to the vertical. 
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1. Introduction 

Disentangling the signatures of various atmospheric waves that 
are supported by structured magnetic atmospheres is a difficult, 
even daunting, task. Recent a dvances in dynamic, 2-D modeling 
of m agnetised atmospheres (Rosenthal et al. 2()02; Bogd an et al.l 
l2003i:,Khomenko & C ollados 2006), highlight the dilemma that 
observers face - i.e., measuring not just singular wave modes 
but the superposition of many magneto-atmospheric modes that 
most likely propagate in different directions and exist in differ- 
ing plasma-y6 environments. This poses a distinct problem, espe- 
cially if studies do not resolve the true size scale of solar mag- 
netic features. 

However, information may be extracted from spatially un- 
resolved structure s by spectropolarim etry (e.g., the penumbral 
flux-tube work of lMiilleretal.1120021) . This approach uses the 
full Stokes polarization spectra (/, Q, U, V), allowing physi- 
cal properties of the emitting plasma to be inferred through the 
application of appropriate model atmospheres. The suitability 
of using Stokes profiles for wave diagnosti cs has been shown 
throtigh numerical sirnulations (se e, e.g., ISolanki & RobertsI 
[l992tMoner & Solanki'l 997", '1999^, while Stokes profiles were 
also used by Rueedi et al. (1998) to interpret magnetic field 
oscillations in a sunspot as resulting from magneto-acoustic- 
gravity waves. 

In this paper we present a method that may identify the form 
of wave which exists in a magnetic environment using the in- 
formation available to full Stokes spectropolarimetry. The ob- 
servational data, their reduction, and details of the form of at- 
mospheric inversion procedure applied are outlined in Sect. |2l 
Results of response function calculations and a Fourier phase 
difference analysis are presented and discussed in Sect. [3] in 
terms of the various forms of wave modes which may exist at 



differing propagation angles, while in Sect.|4]the implications of 
our work are summarized. 



2. Observations 

A small and ellipsoidal sunspot, NOAA 10436, was observed 
on 21 August 2003 using the Tenerife Infrared Polarimeter 
(iMartfnez Fillet et al.|[l999l) attached to the 0.7 m German so- 
lar Vacuum Tower Telescope in Tenerife, Canary Islands. A fast 
scan was performed over the whole sunspot to obtain a global 
picture, comprising 79 slit positions of 0'.'35 width (each with 
a slit integration time of 7 s) incrementally stepped 0'.'4, from 
which the continuum intensity image in Fig. [1] was constructed. 
The sunspot was located somewhat out of disk centre (6 - 26.7°, 
jj. = cos 6 w 0.9) and a light bridge separated the two umbrae. 

Prior to this scan, the slit was fixed across the sunspot and 
the full Stokes vector (/, Q, U, V) was measured (see Fig.|2]for 
example spectra) in a time series comprising 250 stationary-slit 
exposures acquired at a cadence of 14.75 s over 14:39-15:41 UT. 
Seeing conditions were moderate during the observations, with 
an estimated spatial resolution of around r.'5. The most striking 
feature observed in the time series was the oscillatory behaviour 
of Stokes Q, most prominently seen in the inner part of the limb- 
side penumbra (white part of the slit in Fig. [T]!. This oscillation 
in the Stokes Q signal is diplayed in Figs. ^ and|3j/, where a 
dominant 5 min period is observed. 

2. 1 . Spectral Lines 

The recorded spectral region contains two moderately magnet- 
ically sensitive neutral iron lines (Fez 15662.0 A with effective 
Lande factor g = 1.50 and Fei 15665.2 A with g = 0.75) at a 
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Table 1. Atomic parameters of the observed lines. Aq denotes the laboratory wavelength, xi the excitation potential of the lower 
energy level, and log gf the logarithm of the oscillator strength times the multiplicity of the level. The parameters a and cr (in units 
of the Bohr radius, aq) are used to calculate the line broadening by collisions with neutral hydrogen, atoms, while gi, gu, and g are 
the calculated Lande factors of the lower and upper levels, and the effective value, respectively. 
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Fig. 1. Continuum intensity image of NOAA 10436 at 
15666.5±0.5 A. White (black) contours mark the um- 
bral/penumbral (penumbral/quiet Sun) boundaries. The straight 
black line marks the slit position during the time series and the 
white portion indicates the region studied in detail. The dark grey 
pixel within the white portion of the slit marks the position from 
which the profiles in Fig.|2]are taken, while the aiTow in the um- 
bra points to disk centre. 



wavelength sampling of 29.7 mA pixel '. The data reduction 
included polarization demodulation and calibration, flat field- 
ing, dark current and continuum correction, and the removal 
of in strumental cross-talk between the Stokes profiles dCoUadosI 
12003 ). The noise in the reduced data lay below 5 x 10 in units 
of continuum intensity. The wavelength calibration was done 
assuming that the core position of average quiet-Sun intensity 
profiles corresponds to the laboratory wavelength, Aq, shifted 
towards the blue by 500 m s"', the approximate value for the 
blueshif t in these lines deduced from the convective velocity re- 
lation of lNadeaul(ll988h . 

Table [T] presents the atomic data for the spectral lines used 
in this work, where laboratory wavelengths, electro nic config- 
uration s, and excitation potentials were taken from iNave et al] 
d 19941) while the parameters a and cr, which are used in the 
calculation of spectral line broad ening by collisions with n eu- 
tral perturbers, were taken from lAnstee & O'Maral (I1995L sp 
transitions). The two-compone nt model of the quiet Sun from 
iBorrero & Bellot Rubiol (l2002h was used to calcu late empirical 
oscillator strengths for the observed lines, as in iBorrero et alj 
( l2003h . For the 15665 A line, the value for the derived oscil- 
lator strength is especially inaccurate (we estimate +0.2 dex) as 
the intensity profile is blended and the effective quantum num- 
ber for the upper level was beyond the value given in the tables 
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Fig. 2. Example penumbral profiles of Stokes / (a). Stokes Q 
(b). Stokes U (c), and Stokes V (d), each normalized to the local 
continnum intensity, Ic, from the pixel marked dark grey inside 
the region of interest shown in Fig. [T] Wavelengths are relative 
to the 15662 A line, with dotted lines at the blueshift-corrected 
laboratory values. 

of lAnstee & O'Maral (Il995h : here we take the maximum value 
provided by these authors to avoid large extrapolations. 

Figure |2] represents an example of the observed penumbral 
profiles, in this case from the dark grey pixel inside the region 
of interest of Fig. [1] As already mentioned, the intensity pro- 
file of the 15665 A line is blended with an unidentified profile. 
We identify the blend as solar - it varys in strength between 
quiet Sun and umbra - but its origin could not be determined. 
However, it seems that it is not magnetically sensitive because 
even when the field is strong, as in the penumbra, no residual 
polarization signal appears at that wavelength. In our analysis 
we only weakly consider the / spectrum for this line but make 
full use of the polarization spectra (Q, U, V) because, although 
small in magnitude, they provide additional information. 

The circular polarization (Stokes V) shows that the magnetic 
field points downward in the spot but, more interestingly, the lin- 
ear polarization signals (Stokes Q and U) are oppositely signed 
in each line. This is due to the ir opposite Zeeman patterns and 
can be easily proved following iLandi degl'Innocentil (Il992l) . In 
the weak magnetic field limit, it can be written that, 

Q(Ao)^-^-{f + e)Alsm^ycos2xl^] , (1) 
U(Ao)^~{f + e)Alsm^ysm2x[0^_^^ , (2) 
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where y is the inchnation of the magnetic field vector to the Hne- 
of-sight (LOS), X is the azimuthal angle of the magnetic field 
vector in the plane perpendicular to the LOS, Xb is the Zeeman 
wavelength splitting, ^ is the eff'ective Lande factor of the transi- 
tion, and e is a correction factor for anomalous Zeeman triplets 
(7l + ^wJu + \ or + gu) given by. 



e = ^(lUuiJu + D-JLiJL + Df 



-16 [Ju (Ju + 1) + Jl (Jl + 1)] + 4) (gu - gif , 



(3) 



where Jl and Ju are the total angular momentum quantum num- 
bers for the lower and upper transitions, respectively, as defined 
by the electronic configurations (Table [TJ, while and gu are 
the Lande factors of the lower and upper levels whose expres- 
sions can be analyticall y determined as these lines are in pure 
L-S coupling (see, e.g.. Idel Toro Iniestall2003h . Values of ^ H- e 
can be calculated, yielding ^2.2 for Fei 15662 A and ^-1.1 for 
Fe I 15665 A, explaining the opposite linear polarization signals 
observed in Fig. |2] Note, the inversion technique detailed in the 
following section makes use of the full Zeeman splitting pattern 
and does not use the efi'ective Lande factor. 



2.2. Atmospheric Inversion 

The data were inverted using the SPINOR inversion code 
jFru tiger 2000; Frutiger et al. 2000). Prior to the inversion we 
computed the relative Stokes V area asymmetry (iBorrero et alj 
l2004j) . defined as. 



6A 



J V (A) dA 
J \V(A) \dA 



(4) 



Given the small magnitude of 6A observed in the circular polar- 
ization over the inner limb-side penumbra (Fig.[3l7; area between 
the dot-dashed lines) it is likely that these lines are not affected 
by strong gradients in either the magnetic field vector or velocity 
along the L OS. There fore, we have adopted the same model as 
used by Borrero et alj (12004) and referred to as a two-component 
(2-C) inversion. This model consists of two magnetic compo- 
nents and one non-magnetic straylight component - all parame- 
ters except temperature are height independent in each compo- 
nent. The inversion returns the thermodynamic, magnetic, and 
kinematic structure of the atmosphere that provides the best fit 
to the observed Stokes (/, Q, U, V) polarization signals. 

Parameters retrieved from the inner limb-side penumbra 
(Fig- El pixels 1-11) yield a magnetic geometry consisting 
of a near-horizontal component, at local solar inclinations y' ^ 
90 - 125° (55 - 90° from vertical), and a closer-to-vertical one, 
y' « 150 - 16 0° (20 - 30 ° from vertical), consistent with the 
observations of iTitle et al.l (0L993). Note that the retrieved values 
of y', temperature, fi eld strength, a nd filling factor are similar 
to those obtained by iBorrero et alj (l2004i) . From this point on, 
the near-horizontal component will be referred to as flux tube 
(FT) and the less-inclined one as magnetic background (MB) 
follo wing the flux tube inter pretation of Solanki & Montavon 
([1993), Schlichenmaier^sLalJ (fl998l) . iMiiller et all (12002). and 
IBorrero et al. (2005, 20061). Note that this interpreta tion is sub- 
ject to current discussion dSpruit & Scharme"^l2006l) . In princi- 
ple, it may be possible to distinguish between the two scenarios 
by means of an analysis similar to ours, but this requires further 
development of the Spruit and Scharmer scenario and is beyond 
the scope of the current paper. 




Fig. 3. Space-time plots of continuum intensity as a percentage 
of average quiet-Sun continnum (a), magnitude of the relative 
Stokes V area asymmetry (b), and Stokes Q at -1-0.145 A from 
the core of the 15662 A line (c). d): Spatially-averaged Stokes Q 
signal from the region of interest. Only the slit portion extend- 
ing toward solar north east from the northern umbra is shown 
in fl-c. The white (black) contour marks the umbral/penumbral 
(penumbral/quiet Sun) boundary and the dot-dashed lines bound 
the region studied (white area of slit in FiglTJ. 
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Fig. 4. Variation through the limb-side penumbra of the param- 
eters obtained by the inversion at log (r) = : temperature (a), 
field strength (b), local solar inclination (c), and filling factor (d). 
Points mark the temporal mean for each spatial pixel and error 
bars extend over + lcr. 



4 



D. S. Bloomfield et al.: Modified p-modes in penumbral filaments? 






1 












b 


Mog. 


3ockgd. 




8 




F 


ux Tube . 






Mog. Bockgd. fp.b.J 


c 






Flux Tube (p.b.) 




6 








o 










X 


4 








CD 




















D 

cn 


2 






i 


(D 










qI 









i . 




-50 50 100 
Atmospheric Height [km] 



-50 50 100 
Atmospheric Height [km] 



Fig. 5. Height variation of wavelength-integrated 15662 A and 
15665 A Stokes Q velocity response (a), total pressure {b), 
acoustic (c), and Alfvenic {d) wave speeds for the magnetic 
background and flux tube atmospheres (solid and dotted curves, 
respectively). Vertical lines in a mark the COG in each com- 
ponent; vertical dashed and dot-dashed lines in b-d show these 
heights translated into the reference frame of the other atmo- 
sphere by enforcing total pressure balance. 



3. Results & Discussion 

3.1. Height Separations 

In the following analysis the height at which the oscillations are 
observed in both components plays an important role, so that 
we describe in detail how this height is determined. We com- 
pute the line depression response func tion (RF), which is most 
appro priate for the Stokes parameters (iGrossmann-Doerth et al.l 
119881) . The wavelength-integrated RFs of Stokes Q to LOS ve- 
locity are presented in Fig.jSji since the velocity variations were 
most strongly observed in Stokes Q (Fig.O. The RFs for each of 
the components overlap over most of the atmosphere, but their 
center-of-gravity (COG), displayed as vertical lines, reveal a dis- 
tinct separation between the components. These COG heights 
were taken as the origin of the LOS velocity signals. 

One problem with comparing these COG heights is that each 
refers to the height scale of the atmosphere in which the lines 
were computed. The atmospheres returned by the inversion each 
have their own height scale due to the different effective tempera- 
tures (Fig.|4}2) yielding diff'erent pressure scale heights (Fig.lSj?). 
As such, they cannot be simply reduced to a single height scale, 
since pressure balance between the two atmospheres can only be 
enforced at a single height at a time. 

In order to determine vertical height separations between the 
COG heights of the two components this height scale inequality 
must be overcome. This was achieved by enforcing total pressure 
balance between the two atmospheres at the COG heights. For 
example, in Fig.lSj? the value from the FT pressure curve at the 
FT COG (dotted vertical line) is translated onto the MB pressure 
curve, yielding its pressure-balanced COG in the MB reference 
frame (dashed line). Similarly, the value from the MB pressure 
curve at the MB COG (solid vertical line) is translated onto the 
FT pressure curve, providing its pressure-balanced COG in the 
FT reference frame (dot-dashed line). 
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Fig. 6. a): Fourier power spectra from co-spatial magnetic back- 
ground and flux tube atmosphere LOS velocities (solid and dot- 
ted curves, respectively), b): Fourier phase difference spectra be- 
tween the magnetic background and flux tube velocities from the 
eleven analyzed pixels. Darker shading denotes greater Fourier 
coherence and larger symbol size greater cross-spectral power 
c): PDF of phase difference values over the range 2.5 - 4.5 mHz. 
The thick curve displays the measured values and the thin curve 
the best-fit Gaussian profile to the data. 



Through this approach we arrive at representative heights for 
the velocity signals of the MB and FT atmospheres in the refer- 
ence frames of either atmosphere. Values are determined for both 
cases as a consistency check for the method. This allows the ver- 
tical separation distances of the velocity signals to be calculated 
in either of the reference frames - in the MB (FT) atmosphere 
it is the distance between the solid and dashed (dotted and dot- 
dashed) vertical lines. These heights also allow the retrieval of 
characteristic wave propagation speeds from the output inversion 
atmospheres (e.g., Figs.|5}: and|5]i). 



3.2. Time Series Analysis 

All atmospheric parameters except velocity were found to be 
time independent within the error bars. Clear periodic LOS ve- 
locity variations were observed in both atmospheric compo- 
nents, prompting the application of Fourier phase difference 
analysis following Kriigeretal. (2001). The analysis was per- 
formed between the MB and FT velocity signals of eleven pixels 
of the inner limb-side penumbra. The phase difference from one 
pixel is the frequency-dependent phase lag of the FT velocity 
signal with respect to the MB signal, while the coherence mea- 
sures the quality of phase difference variation. 

Figure |6]displays the output from applying such an analysis 
to the observed velocities. Example Fourier power spectra from 
one spatial pixel are given in Fig. ^ where both components 
obviously show power at very similar frequencies. The over- 
plotted phase difference spectra for all eleven of the analyzed 
spatial pixels is presented in Fig. [SJ?, where symbol size rep- 
resents the magnitude of cross-spectral power and shading de- 
notes the coherence. Phase difference values are approximately 
constant in the range showing strongest cross-spectral power 
(~2.5 - 4.5 mHz) and, although close to zero, the probability 
distribution function (PDF) in Fig. |6}: reveals that they are cen- 
tred on -5.5°. This centroid phase difference value was con- 
verted to a time delay between the signals, resulting in values 
of -6.3 s to -3.8 s over the detected range of oscillation fre- 
quencies. Negative phase differences mean that the FT velocity 
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Fig. 7. Variation of propagation angle from the vertical, i^, with 
angular wave mode, I, for acoustic waves at a frequency of 
3.5 mHz. The solid curve shows the case 50 km below the 
equipartition level (where the acoustic and Alfven speeds are 
equal) in a 1.75 kG field inclined at 20°, while the dashed curve 
is 90 km below the equipartition level in a 1.5 kG field inclined 
at 75° (the sampling heights and configurations of the magnetic 
background and flux tube components, respectively). 



leads the MB, agreeing with the COG heights in Fig. |5]for up- 
ward wave propagation. 
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Fig. 8. Comparison of RF-predicted (solid lines) and calculated 
(dotted lines) wave travel distances in the reference frame of the 
magnetic background (left) and the flux tube (right). Cases are 
presented for field-aligned waves propagating at the Alfven (top) 
and slow-mode tube speeds (bottom). The RF-predicted vertical 
height separations have been converted into path length along 
the propagation direction. 



3.3. Wave Modes 

The range of atmospheric parameters existing in a sunspot 
penumbra supports the possibility of many differing forms of 
wave modes. However, all of these could be excited by p-mode 
waves propagating up towards the solar surface and the fact that 
the power peaks at ~3 mHz (Fig.|6]5!) suggests that this is indeed 
the case for the observed wave modes. The dispersion relation 
for acoustic (i.e., p-mode) waves has the form k^c^ - (J- - w^^,, 
where is the sound speed, uj is the angular frequency, and 
is the acoustic cutoff frequency: waves are evanescent for fe^ < 
and propagation occurs for w > Wac- 

Attributing the observed phase differences in the 2.5 - 
4.5 mHz range to propagating waves conflicts with this simple 
picture of evanescent behaviour as the cutoff frequency at the 
photosphere is ~5 mHz in the isothermal case. However, in the 
presence of a magnetic field the acoustic cutoff i s reduced for 
non-v ertical waves in a rather complicated manner & Lerovl 
Il977h . The largest deviation occurs for waves in the strong-field 
limit (i.e., when the Alfven speed, va, is much greater than cs). 
In this situation the cutoff frequency is lowered to Wac cos y' - 
termed the ramp effect - where y' is the magnetic field inclina- 
tion from the vertical. Furthermore, p-modes may travel at an- 
gles away from the vertical at the heights sampled here. This is 
illustrated in Fig. |7] using ray-theory calculations like those of 
ICallvl (I2007h for 3.5 mHz waves in environments similar to our 
two magnetic components: angles of 40 - 60° are possible for 
angular modes I ■x, 760 - 980 around the sampling heights of the 
MB and FT components. 

The situation is further complicated here as these data are not 
recorded in the strong field limit (Fig.|5j > \'a below 70 km 
and 90 km in the MB and FT atm ospheres, respectively; previ- 
ously noted bv lSolanki et alJ[T993l) and the waves are influenced 
by two separate magnetic inclinations. If we assume that inci- 
dent /9-mode waves actually experience some average between 
the differing magnetic inclinations of the MB and FT compo- 
nents then the effective value of field inclination will be in the 



range 40 - 60°. As such, investigation of Fig. 1 in 'Bel & Lerovl 
(Il977i) yields an expected reduction of the cutoff frequency to 
~0.8 a>ac (~4 mHz) for the case where cs * v^, with an absolute 
maximum reduction to 0.5 Wac («2.5 mHz) in the strong field 
limit. However, we note that the concept of a cutoff frequency 
is somewhat ques tionable (see discussions in lSchunker & CalTvi 
l2006l : ICallvll2007l), especially its discrete n ature if radiative cool- 
ing is considered dWebb & Robertslll980h . 

The low-photospheric sampling of the spectral lines means 
that both components are mostly gas dominated over the heights 
sampled (Fig. |5]l. As mentioned previously, differing forms of 
wave can exist, each having certain properties in terms of their 
propagation speed and direction: isotropic acoustic waves prop- 
agate at C5 ; Alfven waves are restricted to the direction of the 
field and move at va; magneto-acoustic slow modes propagate 
along the field at va or, if the sunspot is considered a large "flux 
tube", the tube speed, c\ - c^v^/(c^ -H v^); magneto-acoustic 
fast modes move at, 

2 c\^v\ ^{cl+vlf-4clvl(cosa)^ 
= + 5 , (5) 

where a is the angle between the direction of the field (y') and 
that of wave propagation (tfr). Note that the fast-mode speed dif- 
fers depending on the initial direction of the p-mode waves ex- 
citing these wave modes, the most extreme difference being be- 
tween the case of waves moving toward the sunspot (ingressing; 
a - if/ + y') and those moving away (egressing; a = ip-y') when 
the wave vector occurs in the same azimuthal plane as that of the 
magnetic field. 

Taking these considerations into account, the vertical height 
separations obtained in Sect. 13. H are converted into path lengths 
along the direction of propagation. Probability distribution func- 
tions of RF-predicted path lengths from every space-time pixel 
are presented in Figs. [8] to [12] as solid curves, with values cal- 
culated by combining time delays, wave speeds and propagation 
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Table 2. Differences between COG values of RF-predicted and calculated path length distributions. Values are presented for cases of 
waves excited by p-modes with inclinations, i/r, of 0°, 40°, 50°, and 60° from the vertical. Note that in the case of vertical propagation 
iij/ - 0°) there is no distiction between ingressing or egressing forms of fast-mode waves, hence only one value is provided. When 
sampling a low plasma-y6 environment the Alfven and slow-mode waves are field aligned, remaining invariant to changes in the 
originating p-mode inclination: values are not provided for these two wave modes in the flux-tube reference frame due to their large 
magnitude and obvious incorrectness. 
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Fig. 9. As for Fig. [8] but for acoustic {top) and fast-mode waves Fig. 10. As for Fig. [8] but for acoustic {top) and fast-mode waves 



{bottom) propagating vertically. 



{bottom) propagating at 40° to the vertical. 



angle to the field given as dotted (and in the case of non-vertical 
fast-mode waves also dashed) curves. The difference between 
COG values of the predicted and calculated distributions are 
given in Table|2]as a quantitative measure of the correspondence 
between the various PDFs. 

A comparison between the RF-predicted and calculated path 
lengths of Alfvenic and slow-mode waves is presented in Fig.|8] 
where better correspondence is observed for the case of Alfvenic 
waves over that of slow modes in the MB reference frame. In 
the FT reference frame, however, effectively no correspondence 
is observed between the predicted and calculated PDFs due to 
the large values of field inclination. The PDF comparison for 
these two wave modes does not change for the differing values 
of originating p-mode propagation angle shown in Table |2] as 
the Alfvenic and slow modes remain directed along the mag- 
netic field. The isotropic nature of acoustic and fast-mode waves 
expected in the sampled region of the atmosphere means that al- 
most any angle of propagation could be adopted. An initial con- 
sideration of vertical propagation is provided in Fig. |9] which 
shows that acoustic waves yield a slightly better correspondence 
over fast modes, although neither can be considered successful. 

However, more appropriate values of 40°, 50°, and 60° from 
the vertical are presented based on the previous discussion of 
p-mode behaviour. The case of propagation at 40° to the verti- 
cal is shown in Fig. [10] where it is unclear if acoustic or fast- 



mode waves yield better overlap between predicted and calcu- 
lated PDFs. The distributions given in Fig. [TT| are arrived at 
when considering propagation at 50° to the vertical. In this sce- 
nario fast-mode waves appear to give better correspondence than 
acoustic waves in both the MB and FT reference frames, with the 
case of ingressing fast modes providing an improvement over 
that of egressing fast modes. In the final case studied, for waves 
propagating at 60° to the vertical, all of the PDFs in Fig.[T2]show 
very poor correspondence when compared to the cases of 40° 
and 50° presented above. 

4. Conclusions 

Periodic LOS velocities retrieved from two atmospheric compo- 
nents in a sunspot penumbra have been studied to identify the 
form of wave mode present. The best correspondence between 
RF-predicted and calculated path lengths is observed for fast- 
mode waves propagating toward the sunspot (i.e., ingressing) at 
~50° to the vertical: this scenario has the smallest combination 
of differences between the COG separations of predicted and 
calculated distributions in the two reference frames (Table |2]l. 
The case for fast-mode ingression over egression is supported by 
the location of the sunspot being ~27° from disk centre: egress- 
ing waves in the limb-side (and hence ingressing waves in the 
centre-side) penumbra will have a considerable component of 



D. S. Bloomfield et al.: Modified p-modes in penumbral filaments? 



7 



RF Height Separation 
Time-speed Estimation 



RF Height Seporotion 
Time-speed Estlmotlon 



100 -80 -60 -40 -20 



100 -80 -60 -40 -20 20 



- - Ingressing 
Egressing 




Ingressing 

Egressing 



-100 -80 -60 -40 -20 
Distance [km] 




100 -80 -60 -40 -20 
Distonce [km] 



Fig. 11. As for Fig. [8] but for acoustic (top) and fast-mode waves 
(bottom) propagating at 50° to the vertical. 
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Fig. 12. As for Fig.|8] but for acoustic (top) and fast-mode waves 
(bottom) propagating at 60° to the vertical. 



their plasma motions directed perpendicular to our LOS and thus 
be difficult to observe; ingressing waves in the limb-side (and 
egressing waves in the centre-side) penumbra will be predomi- 
nantly along our LOS. Further support is given by the horizontal 
wavelengths {I x 870, A x 5 Mm or 7") of the p-modes re- 
sponsible for their excitation in the limb-side penumbra: egress- 
ing waves will have traversed two or three "skips" through the 
sunspot, including regions beneath either one or both of the um- 
brae, increasing the likelihood of their absorption or disruption; 
ingressing waves will be on their first "skip" into the outer re- 
gion beneath the spot, remaining relatively strong and coherent. 
Despite the centre-side penumbra lacking strong Stokes Q vari- 
ations like those seen in the limb-side (because of the differing 
magnetic geometries with respect to the LOS), less clearly peri- 
odic velocities of ~5-min period are seen there - compatible with 
ingressing limb-side waves being disrupted after traversing the 
spot. Although these data may not fully validate results from lo- 
cal helioseismology, in regards to wave ingression and egression, 
this may be possible in th e future using 2-0 spec tropolarimetric 
data - e.g., using TESOS dKentischer et all 19981) in VIP mode - 



allowing direct comparison between time-distance analysis and 
this technique. 

This is the first time that spectropolarimetric data have been 
used in this manner to identify a magneto-acoustic wave mode. 
The fact that a fast-mode wave best fits the observational data 
makes qualitative sense as the spectral lines used here sample a 
high-y6 region of the deep photosphere where p-mode waves are 
expected to be modified into fast-mode waves by the presence 
of a magnetic field. As such, this further highlights the role that 
solar internal acoustic waves may play in dynamic phenomena 
both at and above the solar surface. It will be interesting to see 
if the detected form of magneto-acoustic wave changes from the 
case where the velocity response of a spectral line is formed be- 
low the cs - va (i.e., yS » 1) surface to one where it is formed 
above this level. In particular, the ratio of wave amplitude, or en- 
ergy content, observed both above and below the cs - va level 
may help confirm the directi on of the incident p-mode waves 
(c.f.. lSchunker & Callvll200g) . 

This paper illustrates the capability of Stokes spectropo- 
larimetry for improved wave-mode identification over imaging 
studies, which require an assumption about the production of 
intensity variations as well as inferrence of the magnetic field 
geometry that are usually, if at all, provided by potential field ex- 
trapolations from LOS magnetograms. The benefits of the com- 
bined determination of plasma velocities and retrieval of the full 
magnetic vector appear to outweigh the reduction in spatial cov- 
erage caused by using a slit-based instrument. 

Finally, we note that the interpretation may depend to some 
extent on the model employed when carrying out the inver- 
sions. We have restricted ourselves to the simplest such two- 
component model. The use of more sophisticated models could 
lead to further refinements in the results. 
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